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! ! We present preliminary results of the first 7r interferometry (HBT) excitation function 

^ | at intermediate AGS energies. The beam energy evolution of the correlations' dependence 
on rriT, centrality, and emission angle with respect to the reaction plane are discussed. 
Comparisons with predictions of the RQMD cascade model are made. 
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Two-particle intensity interferometry (HBT) measurements have long been used to 



study the geometry and dynamics of heavy ion collisions (see, e.g. [T]]). Pion correlation 
functions are sensitive to the pion source size, shape, decay-time, and long-lived particle 
(e.g. A) production. In addition, dynamic effects such as flow produce space-momentum 
correlations resulting in dependences of the correlation functions on n momentum. 

In this paper, we discuss an excitation function (2-8 AGeV) of n~ HBT measurements. 
Studying the evolution of the correlations as a function of Ef, eam is important for two rea- 
sons. Firstly, a sudden increase, at some Eb eam , in the lifetime of the hadronic fireball has 
long been proposed as a robust signal of the onset of QGP formation . Secondly, the 
sensitivity of correlation functions to the underlying physics makes such measurements 
potent tools to test the dynamics of microscopic models of heavy ion collisions. Many 
models attempt to extrapolate to the RHIC energies. Confidence in the ability to extrap- 
olate (determined by the correct underlying physics and its evolution with energy) would 
be enhanced if the model reproduces an excitation function of detailed HBT systematics. 

Using the large-acceptance EOS Time Projection Chamber |J the E895 collabora- 
tion measured charged particles from Au+Au collisions at 2, 4, 6, and 8 AGeV at the 
Brookhaven AGS. Good particle identification minimized e~ contamination of the n~ 
sample. Momentum resolution, largely due to multiple Coulomb scattering and strag- 
gling in the 3% interaction length target, was on the order of 1.5-3%. The experimental 
correlation functions have been corrected for the momentum resolution with an iterative 
method similar to that employed by the NA44 collaboration M . This correction typically 
increases the fitted A parameter by 15%, and the radii by 5%. 

Track merging and splitting effects were eliminated by a two-track geometrical cut based 
on the tracking algorithm. As expected, this cut discards some pairs (in the "real" and 
event-mixed distributions) at low relative momentum, q. However, due to detector geom- 
etry, this cut preferentially discards low-q pairs at high pr; thus, to minimize phasespace 
bias effects, we restrict our analysis to low pt and use narrow windows in p?. 

A full Coulomb-wave integration over a spatial source of 5-fm Gaussian radius was 
used to generate the Coulomb correction, which was applied pair-wise to the event-mixed 
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denominator of the correlation function. Identical Coulomb correction functions were 
applied to data and to correlation functions from the RQMD (see below). Further details 
of the analysis have been reported previously 0. 

The high quality of the data is seen in Figure [I], where projections in the Bertsch-Pratt 
(BP), or out-side- long, system are shown for midrapidity pions from central events at each 
bombarding energy. The relative momentum q was calculated in the fixed Au+Au cm. 
frame. The functional form 

C(q out , q si de, Qiong) = 1 + Ae-^'^-^ 9 -*-^^ 9 ^- 2 ^ 90 "* 9 ' ^ (1) 

was fit to the data, using a maximum-likelihood technique [IT]. The cross-term [[I Ij R 2 ol 
was consistent with zero in all cases and uncorrelated with the other parameters. 
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Figure 1. Projections of the 3-D 7r~ correlation functions from central collisions 
(~11%<7t). q ou t (top), q S ide (middle), and qi ong (bottom) projections are shown for col- 
lisions at all beam energies measured, y = y cm ±0.35 and pt—0- 1-0.3 GeV/c cuts were 
applied. Projections in a given q-component are integrated over ±30 MeV/c in the other 
components. Superimposed curves show projections of 3-D fits described in text. 

Correlation functions were also constructed in the Yano-Koonin-Podgoretskii (YKP) 
decomposition [1^]; here, the effective lifetime R is fit more directly. We fit to the form 



C(q , q x , m) = 1 + Ae-^-^^-^^^ 



(2) 



The excitation function of the fit results is presented in Figure [2|. Also shown are 
results of fits to correlation functions generated by using the ir~ freeze-out points from 
the RQMD (v2.3) model [13] as input to the two-particle correlator code CRAB [ftl]. Both 
the data and the model show a decrease in the A parameter, due to increased production 
of long-lived 7T~-emitting particles at higher energy 0. 

While the longitudinal radii, Ri ong and R\\, display little dependence on beam energy, 
the observed decrease of the transverse radii R S ide and R out comes as something of a 
surprise. At low energy, the YKP fits suggest that the model produces a tv~ source that is 
too small and too long-lived; this leads in the BP decomposition to an underprediction in 
R S ide, but a reasonable agreement in R ou t, as the space and time effects partially cancel. 
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Part of the reason for the decrease in the apparent transverse size is revealed in Figure |3|. 
It is clear that the m^-dependence of R S ide (and the space-momentum correlation that 
causes it) becomes stronger with E beam in the data, while the model suggests a m T - 
dependence rather independent of energy, and reminiscent of trends at higher energy. 
The observed trend may suggest that collective transverse flow builds with bombarding 
energy, and is only strong enough to affect Rsidei^T) above Eb eam ~4 AGeV. 

Although dynamics determines the m^-dependence of the radii, it is worthwhile to 
check that the HBT radii track somewhat with geometry. Figure [| shows the impact 
parameter (b) dependence of the HBT radii (b was estimated from the charged particle 
multiplicity). A increases with b, as the production of long-lived vr-emitting particles is 
suppressed relative to direct pions. The transverse radii (R S id e and Rout) decrease for 
more peripheral collisions, as expected, while Ri ong shows little 6-dependence. The trends 
suggest that the measured pion source reflects the overlap volume of the colliding nuclei. 

More detailed information may be obtained by studying the HBT signal as a function 
of 7T _ emission angle with respect to the reaction plane. The reaction plane is calculated 
only from momenta of Z< 2 nuclei, for every event |L5|, so auto-correlations are not an 
issue. From an overlap-volume picture, one expects an anisotropic apparent shape in the 
transverse direction for non-central collisions, with a larger spatial scale perpendicular to 
the reaction plane. Deviations may reflect the non-isotropic flow dynamics of the system 



prior to freeze-out, or may carry information concerning the opacity of the source ||16|| . 

Preliminary results for 2, 4, and 6 AGeV collisions at b ~ 5 — 7 fm are shown in Fig. |5|. 
At the lower energy, R S ide (the radius most closely related to geometry [TJ) exhibits a (/in- 
dependence consistent with geometric considerations. While RQMD simulations (with 
perfect reaction plane resolution) display similar trends at all energies, at higher energy, 
the R S ide oscillation is not seen in our data. Since the radii are not corrected for the finite 
dispersion, this is due at least in part to the worsening resolution with which the reaction 
plane is measured. Further study of this novel HBT signal is required. 

In summary, we are mapping out the systematics of pion correlations in the energy 
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Figure 4. Centrality dependence of the BP Figure 5. BP parameters versus tt~ emis- 
fit parameters. sion angle with respect to reaction plane. 



range between the Bevalac and maximum AGS energy. Large jumps in source size or 
lifetime at some collision energy, which might indicate the onset of QGP formation, are 
not observed. Surprisingly, the apparent source size is larger at the lower beam energies; 
this appears largely a consequence of weaker space-momentum correlations there. The 
RQMD model, with or without meanfield effects, does not reproduce the data; in the 
model at low energy, the effective size is too small, the lifetime too large, and R(mr) does 
not evolve with beam energy. The impact parameter dependence of the radii follows naive 
expectations from geometry. An HBT analysis correlated with the event-wise reaction 
plane reveals a significant oscillation in R S id e , at low beam energy. 
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